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Abstract 

Vortex line and magnetic line representations are introduced for description of 
flows in ideal hydrodynamics and MHD, respectively. For incompressible fluids it 
is shown that the equations of motion for vorticity fi and magnetic field with the 
help of this transformation follow from the variational principle. By means of this 
representation it is possible to integrate the system of hydrodynamic type with 
the Hamiltonian = / |ri|dr. It is also demonstrated that these representations 
allow to remove from the noncanonical Poisson brackets, defined on the space of 
divergence-free vector fields, degeneracy connected with the vorticity frozenness for 
the Euler equation and with magnetic field frozenness for ideal MHD. For MHD a 
new Weber type transformation is found. It is shown how this transformation can 
be obtained from the two-fluid model when electrons and ions can be considered as 
two independent fluids. The Weber type transformation for ideal MHD gives the 
whole Lagrangian vector invariant. When this invariant is absent this transforma- 
tion coincides with the Clebsch representation analog introduced by Zakharov and 
Kuznetsov. 
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1 Introduction 



There are a large number of works devoted to the Hamiltonian description of the ideal 
hydrodynamics (see, for instance, the review and the references therein). This ques- 
tion was first studied by Clebsch (a citation can be found in Ref. P]), who introduced 
for nonpotential flows of incompressible fluids a pair of variables A and /i (which later 
were called as the Clebsch variables). A fluid dynamics in these variables is such that 
vortex lines represent themselves intersection of surfaces A = const and fi = const and 
these quantities, being canonical conjugated variables, remain constant by fluid advection. 
However, these variables, as known (see, i.e.,0]) describe only partial type of flows. If A 
and fi are single-valued functions of coordinates then the linking degree of vortex lines 
characterizing by the Hopf invariant occurs to be equal to zero. For arbitrary flows 
the Hamiltonian formulation of the equation for incompressible ideal hydrodynamics was 
given by V.I. Arnold 0. The Euler equations for the velocity curl O = curlv 



^— = curl V X III 

dt ^ ^ 



divv = 



are written in the Hamiltonian form, 



'dt 



{n,n}, 



by means of the noncanonical Poisson brackets 



{F,G} 




'5F , 5G 
curl X curl — 
oil oil 



dr 



where the Hamiltonian 



XI) 



:i.2) 



;i.3) 



;i-4) 



coincides with the total fluid energy. 

In spite of the fact that the bracket (|1.3|) allows to describe flows with arbitrary 
topology its main lack is a degeneracy. By this reason it is impossible to formulate the 
variational principle on the whole space S of divergence-free vector fields. 

The cause of the degeneracy, namely, presence of Casimirs annulling the Poisson 
bracket, is connected with existence of the special symmetry formed the whole group 
- the relabeling group of Lagrangian markers (for details see the reviews [|, 0). All 
known theorems about the vorticity conservation (the Ertel's, Cauchy's and Kelvin's the- 
orems, the frozenness of vorticity and conservation of the topological Hopf invariant) are 
a sequence of this symmetry. The main of them is the frozenness of vortex lines into fluid. 
This is related to the local Lagrangian invariant - the Cauchy invariant. The physical 
meaning of this invariant consists in that any fluid particle remains all the time on its 
own vortex line. 

The similar situation takes place also for ideal magneto-hydrodynamics (MHD) for 
barotropic fluids: 

Pt + V(pv) = 0, (1.5) 



Vi + (vV)v = —Vw{p) + - — [curlh x h] 

47rp 



;i.6) 
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hj = curl[v X h]. (1.7) 

Here p is a plasma density, w{p) plasma entalpy, v and h are velocity and magnetic fields, 
respectively. As well known (see, for instance, |P-|T^), the MHD equations possesses 
one important feature - frozenness of magnetic field into plasma which is destroyed only 
due to dissipation (by finite conductivity). For ideal MHD combination of the continuity 
equation (|1.5| ) and the induction equation (|1.7| ) gives the analog of the Cauchy invariant 
for MHD. 



The MHD equations of motion ( |1.5[ - |1.7| ) can be also represented in the Hamiltonian 
form, 

Pt = {p,n} ht = {h,n}, v, = {v,7^}, 

by means of the noncanonical Poisson brackets |T^ 



(1.8) 




6G 

curl—— X -— 
on ov 



6F 

curl— X — 
on ov 



;i.9) 



6£ 6G 
(5v (5v 




This bracket is also degenerated. For instance, the integral / (v, h)dr, which characterizes 
mutual linkage knottiness of vortex and magnetic lines, is one of the Casimirs for this 
bracket. 

The analog of the Clebsch representation in MHD serves a change of variables sug- 
gested in 1970 by Zakharov and Kuznetsov |]l|: 



V0 + 



[h X curl S] 
P 



;i.io) 



New variables (0, p) and h, S represent two pairs canonically conjugated quantities with 
the Hamiltonian coinciding with the total energy 



n 



V 

T 



p— + pe{p) + — \ dr. 



In the present paper we suggest a new approach of the degeneracy resolution of the 
noncanonical Poisson brackets by introducing new variables, i.e., Lagrangian markers 
labeling each vortex lines for ideal hydrodynamics or magnetic lines in the MHD case. 

The basis of this approach is the integral representation for the corresponding frozen- 
in field, namely, the velocity curl for the Euler equation and magnetic field for MHD. We 
introduce new objects, i.e., the vortex lines or magnetic lines and obtain the equations 
of motion for them. This description is a mixed Lagrangian-Eulerian description, when 
each vortex (or magnetic) line is enumerated by Lagrangian marker, but motion along 
the line is described in terms of the Eulerian variables. Such representation removes 
all degeneracy from the Poisson brackets connected with the frozenness, remaining the 
equations of motion to be gauge invariant with respect to re-parametrization of each line. 
Important, that the equations for line motion, as the equations for curve deformation, are 
transverse to the line tangent. 

It is interesting that the line representation also solves another problem - the equations 
of line motion follow from the variational principle, being Hamiltonian. 
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This approach allows also simply enough to consider the limit of narrow vortex (or 
magnetic) lines. For two-dimensional flows in hydrodynamics this "new" description 
corresponds to the well-known fact, namely, to the canonical conjugation of x and y 
coordinates of vortices (see, for instance, 0). 

The Hamiltonian structure introduced makes it possible to integrate the three-dimen- 
sional Euler equation ( |1.2| ) with Hamiltonian Ti. = J \ n\dr. In terms of the vortex lines 
the given Hamiltonian is decomposed into a set of Hamiltonians of noninteracting vortex 
lines. The dynamics of each vortex lines is, in turn, described by the equation of a vortex 
induction which can be reduced by the Hasimoto transformation [jl5| to the integrable 
one-dimensional nonlinear Schrodinger equation. 

For ideal MHD a new representation - analog of the Weber transformation - is found. 
This representation contains the whole vector Lagrangian invariant. In the case of ideal 
hydrodynamics this invariant provides conservation of the Cauchy invariant and, as a 
sequence, all known conservation laws for vorticity (for details see the review 0). It 
is important that all these conservation laws can be expressed in terms of observable 
variables. Unlike the Euler equation, these vector Lagrangian invariants for the MHD 
case can not be expressed in terms of density, velocity and magnetic field. It is necessary 
to tell that the analog of the Weber transformation for MHD includes the change of 
variables (|1.10| ) as a partial case. The presence of these Lagrangian invariants in the 
transform provides topologically nontrivial MHD flows. 

The Weber transform and its analog for MHD play a key role in constructing the vortex 
line (or magnetic line) representation. This representation is based on the property of 
frozenness. Just therefore by means of such transform the noncanonical Poisson brackets 
become non-degenerated in these variables and, as a result, the variational principle may 
be formulated. Another peculiarity of this representation is its locality, establishing the 
correspondence between vortex (or magnetic) line and vorticity (or magnetic field). This 
is a specific mapping, mixed Lagrangian-Eulerian, for which Jacobian of the mapping can 
not be equal to unity for incompressible fluids as it is for pure Lagrangian description. 



2 General remarks 

We start our consideration from some well known facts, namely, from the Lagrangian 
description of the ideal hydrodynamics. 

In the Eulerian description for barotropic fluids, pressure p = p{p), we have coupled 
equations - discontinuity equation for density p and the Euler equation for velocity: 

Pi + divpv = 0, (2.1) 

vt + (vV)v = -Vw{p), dw{p) = dp/ p. (2.2) 

In the Lagrangian description each fluid particle has its own label. This is three-dimen- 
sional vector a, so that particle position at time t is given by the function 

x = x(a,t). (2.3) 

Usually initial position of particle serves the Lagrangian marker: a = x(a,0). 
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In the Lagrangian description the Euler equation ( p.2|) is nothing more than the New- 
ton equation: 

X = — Vw. 

In this equation the second derivative with respect to time t is taken for fixed a, but the 
r.h.s. of the equation is a function of t and x. Exchiding from the latter the x-dependence, 
the Euler equation takes the form: 

= (2.4) 

oak oak 

where now all quantities are functions of t and a. 

In the Lagrangian description the continuity equation (|2.1| ) is easily integrated and 
the density is given through the Jacobian of the mapping (|2.3| ) J = det{dxi/dak): 



P = (2.5) 

Now let us introduce a new vector, 

Uk = T^Vi , (2.6) 
oak 

which has a meaning of velocity in a new curvilinear system of coordinates or it is possible 
to say that this formula defines the transformation law for velocity components. It is worth 
noting that (|2.6| ) gives the transform for the velocity v as a co-vector. 

The straightforward calculation gives that the vector u satisfies the equation 

duk 9 /v2 \ 

-dr = d^k[Y-'')- ^2-^) 

In this equation the right-hand-side represents gradient relative to a and therefore the 
"transverse" part of the vector u will conserve in time. And this gives the Cauchy invari- 
ant: 

-curl„u = 0, (2.8) 

at 

or 

curUu = I. (2.9) 

If Lagrangian markers a are initial positions of fiuid particles then the Cauchy invari- 
ant coincides with the initial vorticity: I = rio(a). This invariant is expressed through 
instantaneous value of fi(x, t) by the relation 

no(a) =J(ri(x,t)V)a(x,t) (2.10) 



where a = a(x,t) is inverse mapping to (|2.3| ). Following from ( |2.1CI|) relation for B = fi/p, 

Oa ■ 

Boi{a) = --!-Bk{x,t), 

OXk 

shows that, unlike velocity, B transforms as a vector. 
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By integrating the equation ( p.7|) over time t we arrive at the so-called Weber trans- 
formation 

u(a,t) = uo(a)+Va$, (2.11) 
where the potential $ obeys the Bernoulli equation: 

^ = y--(p) (2.12) 

with the initial condition: ^\t=o = 0. For such choice of $ a new function uo(a) is 
connected with the "transverse" part of u by the evident relation 

curia uo (a) = I. 

The Cauchy invariant I characterizes the vorticity frozenness into fluid. It can be got 
by standard way considering two equations - the equation for the quantity B = fl/p, 

c/B .„ , 

- = (BV)v. (2.13) 

and the equation for the vector 5x = x(a + Sa) — x(a) between two adjacent fluid parti- 
cles: 

^ = (5xV)v, (2.14) 

The comparison of these two equations shows that if initially the vectors 5x are parallel 
to the vector B, then they will be parallel to each other all time. This is nothing more 
than the statement of the vorticity frozenness into fluid. Each fluid particle remains all 
the time at its own vortex line. The combination of Eqs. (|2.13|) and (p.l4|) leads to the 



Cauchy invariant. To establish this fact it is enough to write down the equation for the 
Jacoby matrix Jij = dxi/daj which directly follows from (|2.14| ): 



one 



d dtti dtti dvj 
dt dxk dxj dxk ' 

that in combination with Eq. (|2.13| ) gives conservation of the Cauchy invariant ( p.9|) . 

If now one comes back to the velocity fleld v then by use of Eqs. ( p.6|) and ( p. 11 
can get that 

V = MofcVafc + (2.15) 

where gradient is taken with respect to x. Here the equation for potential $ has the 
standard form of the Bernoulli equation: 

$i + (vV)<I>-y +W(p) =0. 

It is interesting to note that relations (|2.9| ), as equations for determination of x(a,t), 
unlike Eqs ( |2.7| ), are of the flrst order with respect to time derivative. This fact also 
reflects in the expression for velocity (|2.15|) which can be considered as a result of the 



partial integration of the equations of motion (|2.7|) . Of course, the velocity fleld given 



by ( |2.15| ) contain two unknown functions: one is the whole vector a(x, t) and another 
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is the potential $. For incompressible fluids the latter is determined from the condition 
divv = 0. In this case the Bernoulli equation serves for determination of the pressure. 

Another important moment connected with the Cauchy invariant is that it follows 
from the variational principle (written in terms of Lagrangian variables) as a sequence of 
relabelling symmetry remaining invariant the action (for details, see the reviews 0,0). 
Passing from Lagrangian to Hamiltonian in this description we have no any problems with 
the Poisson bracket. It is given by standard way and does not contain any degeneracy 



that the noncanonical Poisson brackets (|1.3|) and ( |1.9|) have. One of the main purposes of 
this paper is to construct such new description of the Euler equation (as well as the ideal 
MHD) which, from one side, would allow to retain the Eulerian description, as maximally 
as possible, but, from another side, would exclude from the very beginning all remains 
from the gauge invariance of the complete Euler description connected with the relabeling 
symmetry. 

As for MHD, this system in one point has some common feature with the Euler 
equation: it also possesses the frozenness property. The equation for h/p coincides with 
( p.l3| ) and therefore dynamics of magnetic lines is very familiar to that for vortex lines of 
the Euler equation. However, this analogy cannot be continued so far because the equation 
of motion for velocity differs from the Euler equation by the presence of pondermotive 
force. This difference remains also for incompressible case. 



3 Vortex line representation 

Consider the Hamiltonian dynamics of the divergence-free vector field ft{r,t), given by 
the Poisson bracket (|1.3|) with some Hamiltonian Ti. Q: 



— — = curl 
dt 



curl — - X il 



(3.1) 



As we have said, the bracket ( |1.3| ) is degenerate, as a result of which it is impossible 
to formulate the variational principle on the entire space S of solenoidal vector fields. 
It is known that Casimirs /, annulling Poisson brackets, distinguish in S invariant 
manifolds Ai / (symplectic leaves) on each of which it is possible to introduce the standard 
Hamiltonian mechanics and accordingly to write down a variational principle. We shall 



show that solution of this problem for the equations (|3.1| ) is possible on the base of the 
property of frozenness of the field fl{r, t), which allows to resolve all constrains, stipulated 
by the Casimirs, and gives the necessary formulation of the variational principle. 

To each Hamiltonian Ti. - functional of 0(r,t) - we associate the generalized velocity 

v(r)=curl^. (3.2) 

However one should note that the generalized v(r) is defined up to addition of the vector 
parallel to O: 

•^1^^^ T7T ~* curl — + all, 
oil oil 

^The Hamiltonian ( |1.4|) corresponds to ideal incompressible hydrodynamics. 
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that in no way does change the equation for fl. Under the condition [ft ■ Va) = a new 
generahzed velocity will have zero divergence and the frozenness equation can be 
written already for the new v(r). A gauge changing of the generalized velocity corresponds 
to some addition of a Casimir to the Hamiltonian : 

n^n + f; {/,..} = 0. 

Hence becomes clear that the transformation 

X = x(a, t) 

of the initial positions of fluid particles x(a, 0) = a by the generalized velocity field v(r) 
through solution of the equation 

X = v(x,t) (3.3) 

is defined ambiguously due to the ambiguous definition of v(r) by means of ( [3.21 ). There- 
fore using full Lagrangian description to the systems (|3.1|) becomes ineffective. 

Now we introduce the following general expression for fi(r), which is gauge invariant 
and fixes all topological properties of the system that are determined by the initial field 

n{r, t) = J 5{r - R(a, t))(^^o(a) Va)R(a, t)d^a. (3.4) 

Here now 

r = R(a,t) (3.5) 

does not satisfy any more the equation (|373|) and, consequently, the mapping Jacobian 
J = det||9R/(9a|| is not assumed to equal 1, as it was for full Lagrangian description of 
incompressible fluids. 

It is easily to check that from condition (Vaf^o(a)) = it follows that divergence of 
( ^.4| ) is identically equal to zero. 
The gauge transformation 

R(a) ^R(an„(a)) (3.6) 

leaves this integral unchanged if a^^g is arisen from a by means of arbitrary nonuniform 
translations along the field line of r2o(a). Therefore the invariant manifold Ain^ of the 
space S, on which the variational principle holds, is obtained from the space 7^ : a ^ R of 
arbitrary continuous one-to-one three-dimensional mappings identifying TZ elements that 
are obtained from one another with the help of the gauge transformation ( p.6| ) with a 
fixed solenoidal field r2o(a). 

The integral representation for f2 ( p.4| ) is another formulation of the frozenness con- 
dition - after integration of the relation (|3.4| ) over area a, transverse to the lines of fl, 
follows that the flux of this vector remains constant in time: 



Ja(t) Ja(0) 



la(t) Ja(0) 

Here a{t) is the image of a{0) under the transformation ( p.5| ). 

It is important also that rio(a) can be expressed explicitly in terms of the instantaneous 
value of the vorticity and the mapping a = a(r,t), inverse to (|3.5|) . By integrating over 
the variables a in the relation 
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where r2o(a) can be represented in the form: 



ftoia) = det||9R/(9a||(rj(r)V)a. (3.8) 

This formula is nothing more than the Cauchy invariant ( p.9|) . We note that according to 
Eq. ( |3.7| ) the vector b = (r2g(a) Va)R(a) is tangent to f2(R). It is natural to introduce 
parameter s as an arc length of the initial vortex lines f2o(a) so that 

In this expression depends on the transverse parameter u labeling each vortex line. In 
accordance with this, the representation ( p.4|) can be written in the form 

rj(r, t)= J ^Q{u)(fu J 6{r - R(s, z/, t))^ds, (3.9) 

whence the meaning of the new variables becomes clearer: To each vortex line with index 
z/ there is associated the closed curve 

r = R(s, z/, t), 

and the integral (|3.9| ) itself is a sum over vortex lines. We notice that the parametrization 
by introduction of s and u is local. Therefore as global the representation ( |3.9| ) can be 
used only for distributions with closed vortex lines. 

To get the equation of motion for R(z/,s, t) the representation ( PTP| ) (in the general 
case - ( p. 4]) ) must be substituted in the Euler equation (|3.1| ) and then a Fourier transform 



with respect to spatial coordinates performed. As a result of simple integration one can 
obtain: 

k X y no{iy)d^u J dse-'^^[Rs x {Rt(z/,s,t) - v(R,t)}] = 0. 

This equation can be resolved by putting integrand equal identically to zero: 

[R, X Kt{u, s, t)] = [R, X v(R, t)]. (3.10) 

With this choice there remains the freedom in both changing the parameter s and re- 
labelling the transverse coordinates u. In the general case of arbitrary topology of 
the field r2o(a) the vector R^ in the equation (|3.1CI| ) must be replaced by the vector 



b = (r2o(ci)Va)R(a, t). Notice that, as it follows from ( p.lOp and (|3.7|), a motion of a 
point on the manifold Aiuo is determined only by the transverse to fi(r) component of 
the generalized velocity. 

The obtained equation ( p.lOj ) is the equation of motion for vortex lines. In accordance 
with ( ^.1U| ) the evolution of each vector R is principally transverse to the vortex line. The 
longitudinal component of velocity does not effect on the line dynamics. 

The description of vortex lines with the help of equations ( p. 9] ) and ( ^.lU ) is a mixed 



Lagrangian-Eulerian one: The parameter u has a clear Lagrangian origin whereas the 
coordinate s remains Eulerian. 
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4 Variational principle 



The key observation for formulation of the variational principle is that the following 
general equality holds for functionals that depend only on O: 



b X curl 



6F 



sn(R) 



6F 



(5R(a) 



On 



(4.1) 



For this reason, the right-hand-side of ( p.lOp equals the variational derivative STi/SH: 

6n{n{R}} 



[(rio(a)Va)R(a) x Ri(a)] 



(5R(a) 



no 



(4.2) 



It is not difficult to check now that the equation ([4.2| ) described dynamics of vortex line 
is equivalent to the requirement of extremum of the action [6S = 0) with the Lagrangian 

C = - y d='a(|R,(a) X R(a)| • (n„(a)VJR(a)) - H({n{R}}). (4.3) 

Thus, we have introduced the variational principle for the Hamiltonian dynamics of 
the divergence-free vector field topologically equivalent to fio(ci)- 

Let us discuss some properties of the equations of motion ( [4.2|) , which are associated 
with excess parametrization of elements of A^^o by objects from TZ. We want to pay 
attention to the fact that From Eq. ( [4.1|) follows the property that the vector b and 
5F/5R(a) are orthogonal for all functionals defined on Ain^. In other words the varia- 
tional derivative of the gauge-invariant functionals should be understood (specifically, in 
(O) as 

6F 



P 



5R(a) ' 



where P,, 



is a projector and r = b/|b| a unit tangent (to vortex line) 



vector. Using this property as well as the transformation formula (|4.1| ) it is possible, 
by a direct calculation of the bracket (|1.3|) , to obtain the Poisson bracket (between two 
gauge-invariant functionals) expressed in terms of vortex lines: 



b- 



p 



6F 



X P- 



6G 



5R(a) 5R(a 



(4.4) 



The new bracket ( |4.4| ) does not contain variational derivatives with respect to r2o(a). 
Therefore, with respect to the initial bracket the Cauchy invariant r2o(a) is a Casimir 
fixing the invariant manifolds Aiuo which it is possible to introduce the variational 
principle ( [4.3|) . 

In the case of the hydrodynamics of a superfiuid liquid a Lagrangian of the form 
( [4.3|) was apparently first used by Rasetti and Regge [0 to derive an equation of motion, 
identical to Eq. ( ^.l(j| ), but for a separate vortex filament. Later, on the base of the results 
17|, Volovik and Dotsenko Jr. [18| obtained the Poisson bracket between the coordinates 
of the vortices and the velocity components for a continuous distribution of vortices. The 
expression for these brackets can be extracted without difficulty from the general form 
for the Poisson brackets (J4.4|) . However, the noncanonical Poisson brackets obtained in 
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T7| , |18] must be used with care. Their direct apphcation gives for the equation of motion 



of the coordinate of a vortex filament an answer that is not gauge-invariant. For a general 
variation, which depends on time, additional terms describing flow along a vortex appear 
in the equation of motion. For this reason, the dynamics of curves (including vortex lines) 
is in principle "transverse" with respect to the curve itself. 

We note that for two-dimensional (in the x — y plane) flows the variational principle 
for action with the Lagrangian ( |4.3| ) leads to the well-known fact that X(z/, t)- and Ylu, t)- 
coordinates of each vortex are canonically conjugated quantities (see M). 



5 Integrable hydrodynamics 



Now we present an example of the equations of the hydrodynamic type ( |3.1j ), for which 
transition to the representation of vortex lines permits to establish of the fact of their 
integrability |T^. 

Consider the Hamiltonian 

n{n{r)} = [ \n\dr (5.1) 



and the corresponding equation of frozenness ( p.l|) with the generalized velocity 

V = curl (n/n) . 



We assume that vortex lines are closed and apply the representation ( |3.9| ). Then due to 
( p.7|) the Hamiltonian in terms of vortex lines is decomposed as a sum of Hamiltonians of 
vortex lines: 

n{R} = I \n^{v)\d^v ^ 



ds 



ds. 



(5.2) 



The standing here integral over s is the total length of a vortex line with index v. Accord- 
ing to (|4.2| ), with respect to these variables the equation of motion for the vector R(z/, s) 
is local, it does not contain terms describing interaction with other vortices: 



?7[r X Rt(i/,s,t)] = [r X [r X 



(5.3) 



Here rj = sign(fio), t = R-s/lH-sl is the unit vector tangent to the vortex line. 

This equation is invariant against changes s s{s,t). Therefore the equation 
can be resolved relative to Rt up to a shift along the vortex line - the transformation 
unchanged the vorticity ft. This means that to find O it is enough to have one solution 
of the equation 

?7|R,|Rt= [rxrJ+/?R„ (5.4) 



which follows from ( ^.3|) for some value of /?. Arisen from here equation for r as a function 
of filament length / {dl = \Ils\ds) and time t (by choosing a new value /3 = 0) reduces to 
the integrable one- dimensional Landau- Lif shits equation for a Heisenberg ferromagnet: 



dr 



T X 
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This equation is gauge-equivalent to the ID nonhnear Schrodinger equation and, 
for instance, can be reduced to the NLSE by means of the Hasimoto transformation p!5|: 



where K{l,t) is a curvature and xi^^t) the hne torsion. 

The considered system with the Hamiltonian (|5.1|) has direct relation to hydrodynam- 



ics. As known (see the paper [15| and references therein), the local approximation for thin 
vortex filament (under assumption of smallness of the filament width to the characteristic 
longitudinal scale) leads to the Hamiltonian (|5.2| ) but only for one separate line. Respec- 
tively, the equation (|3.1|) with the Hamiltonian (|5.1| ) can be used for description of motion 
of a few number of vortex filaments, thickness of which is small compared with a distance 
between them. In this case (nonlinear) dynamics of each filament is independent upon 
neighbor behavior. In the framework of this model singularity appearance ( intersection 
of vortices) is of an inertial character very similar to the wave breaking in gas-dynamics. 
Of course, this approximation does not work on distances between filaments comparable 
with filament thickness. 

It should be noted also that for the given approximation the Hamiltonian of vortex 
line is proportional to the filament line whence its conservation follows that, however, 
in no cases is adequate to behavior of vortex filaments in turbulent flows where usually 
process of vortex filament stretching takes place. It is desirable to have the better model 
free from this lack. A new model must necessarily describe nonlocal effects. 

In addition we would like to say that the list of equations ( |3.1| ) which can be integrated 
with the help of representation (|3.9|) is not exhausted by ( ^.1|) . So, the system with the 
Hamiltonian 

n^Wr)} = I \n\xdv (5.5) 

is gauge equivalent to the modified KdV equation 

3 

^I't + ^iii + = - 

the second one after NLSE in the hierarchy generated by Zakharov-Shabat operator. As 
against previous model (|5.1|) some physical application of (|5.5|) has not yet been found. 



6 Lagrangian description of MHD 

Consider now how the relabelling symmetry works in the ideal MHD. First, rewrite equa- 
tions of motion ( |1.5[ - |1.7D in the Lagrangian representation by introducing markers a for 



fiuid particles 

X = x(a, t) 

with 

v(x, t) = x(a, t). 

In this case the continuity equation ( [1.5[ ) and the equation for magnetic field ( |1.7| ) can 
be integrated. The density and the magnetic field are expressed in terms of the Jacoby 
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matrix by means of Eq. ( p.5|) and by the equation 

dx. 



Bi{x,t) 



dai 



(6.1) 



where B = h/p. In the latter transformation the Jacoby matrix serves the evolution op- 
erator for vector B. The vector B, in turn, transforms as a vector. 

In terms of Lagrangian variables the equation of motion (|1.6| ) is written as follows 



dx 



Xi 



dw{p) 



+ 



J dx ■ 

■ [curl h X h] 



With the help of relation ( |6.1|) and Eq. 
equation 



the vector u given by ( p.6|) will satisfy the 



du 

~dt 



~w \ - — [Bo (a) X curl„H] 



(6.3) 



Here vector Bo(a) = ho(a) / po{3i) is a Lagrangian invariant and H represents the co-adjoint 
transformation of the magnetic field, analogous to 



dx„ 
(9a,- 



Now by analogy with ( pj.7| ) and ( ^.llj ), integration of Eq.( |6.3| ) over time leads to the Weber 
type transformation: 



u(a,t) = uo(a) + Va$ + Bo(a) x curUS 



(6.4) 



Here uo(a) is a new Lagrangian invariant which can be chosen as pure transverse, namely, 
with diVa uq = 0. This new Lagrangian invariant cannot be expressed through the observed 
physical quantities such as magnetic field, velocity and density. In spite of this fact, as 
it will be shown in the next section, the vector Lagrangian invariant uo(a) has a clear 
physical meaning. As far as new variables $ and S, they obey the equations: 



d<l> 
lit 



V 

T 



dS H _ , 

The transformation (O) for velocity v(x, t) takes the form: 

"h 



X curl S 



P 



V =Mofc(a)Vafc + V$ + 

where S is the vector S transformed by means of the rule 

Si{x,t) 



(6.5) 



dak ^ I ,^ 



dx 

In Eulerian description $ satisfies the Bernoulli equation 



w = 



(6.6) 
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and equation of motion for S is of the form: 

dS h 

^ + 1 [v X curlS] + Vtpi = 0. (6.7) 

ot 47r 

For uq = the transformation (|6.5| ) was introduced for ideal MHD by Zakharov and 
Kuznetsov in 1970 In this case magnetic field h and vector S as well as $ and p 
are two pairs of canonically conjugated variables. It is interesting to note that in the 
canonical case the equations of motion for S and $ obtained in |T]] coincide with ( |6.6| ) 
and ( |6.7|) . However, the canonical parametrization describes partial type of flows, in 
particular, it does not describe topological nontrivial flows for which mutual knottiness 
between magnetic and vortex lines is not equal to zero. This topological characteristics 
is given by the integral /(v, h)(ix. Only when uq 7^ this integral takes non-zero values. 



7 Frozen-in MHD fields 

To clarify meaning of new Lagrangian invariant uo(a) we remind that the MHD equations 



(1.5-L7D can be obtained from two-fluid system where electrons and ions are considered 
as two separate fluids interacting each other by means of self-consistent electromagnetic 
field. The MHD equations follow from two-fluid equations in the low-frequency limit 
when characteristic frequencies are less than ion gyro-frequency. The latter assumes i) 
neglecting by electron inertia, ii) smallness of electric field with respect to magnetic field, 
and iii) charge quasi-neutrality. We write down at first some intermediate system called 
often as MHD with dispersion 



Attc 

curl curlA = {uiVi — n2V2), (7.1) 

c 

{dt + viV)mvi = - {-At + [vi X curl A]) - V—, (7.2) 

p 3f 
= — (-At + [v2xcurl A])-V— . (7.3) 

c on2 

In these equations A is the vector potential so that the magnetic field h = curlA and 
electric field E = —-At. This system is closed by two continuity equations for ion density 
ni and electron density 

ni^t + V(nivi) = 0, n2,t + V(n2V2) = 0. (7.4) 

In this system vi 2 are velocities of ion and electron fluids, respectively. The first equation 
of this system is a Maxwell equation for magnetic field in static limit. The second equation 
is equation of motion for ions. The next one is equation of motion for electrons in which we 
neglect by electron inertia. By means of the latter equation one can obtain the equation 
of frozenness of magnetic field into electron fluid (this is another Maxwell equation): 

= curl[v2 X h]. 
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Applying the operator div to (|7.1|) gives with account of continuity equations the quasi- 
neutrahty condition: Ui = n2 = n. Next, by excluding n2 and V2 we have finally the 
MHD equations with dispersion in its standard form pO| : 

1 d 

idt + vV)mv = —Vwin) H [curl h x h], win) = —e(n,n), 

Ann on 

nt + V(nv) = 0, h( = curl ( v — curl h ) x h , (7.5) 

LV Anen J J 

where vi = v, and e{n,n) is internal energy density so that w{n) is entalpy per one 

pair ion-electron. The classical MHD follows from this system in the limit when the last 

term c/(47ren)curlh in equation ( [7.5|) should be neglected with respect to v. At the same 

time, the vector potential A must be larger characteristic values of (mc/e)v in order to 

provide inertia and magnetic terms in Eq. ( |7.2| ) being of the same order of magnitude. 

Both requirements are satisfied if e = c/{ujpiL) << 1 where L is a characteristic scale of 



magnetic field variation and cUpi = yAirne'^/m is ion plasma frequency. 

Unlike MHD equations dT^jT^), the given system has two frozen-in fields. These are 
the field fio = — frozen into electron fiuid and the field 



frozen into ion component: 



where 



Cti = curl(v + —A) = - O2 

mc 



Qit = curl [v X ill] 

= curl [v2Xf22] 



V2 = V — curl h. 

Anen 



Hence for both fields one can construct two Cauchy invariants by the same rule (|2.9| ) as 
for ideal hydrodynamics: 

Oio(a) = Ji(f^i(x, t)V)a(x,t) (7.6) 
where a(x,)f:) is inverse mapping to x = Xi(a, t) which is solution of the equation x= v(x,t); 

f^2o(a2) = J2(^^2(x,t)V)a2(x,t) (7.7) 

with a2(x,t) inverse to the mapping x = X2(a2,t) and x = V2(x,t). 

In order to get the corresponding Weber transformation for MHD as a limit of the 
system it is necessary to introduce two momenta for ion and electron fiuids: 

= rnv+-A (7.8) 
c 

P2 = --A. (7.9) 

c 

In these expressions the terms containing the vector potential are greater sum of pi and p2 
in parameter e. For each momentum in Lagrangian representation one can get equations, 
analogous to (U), (^) : 

dxk dpik 



dan dt 

dxk dp2k 
da2i dt 




(7.10) 
(7.11) 
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By introducing the vector p for each type of fluids, by the same rule as (|2.6| 



Vi 



dxk 
da. 



after integration over time of equations of motion for p one can arrive at two Weber 
transformations for each momentum: 



Pi = Pii{ai)Vaii + V^i, 

P2 = P2i(a2)Va2i + V$2- 



(7.12) 
(7.13) 



In the limit e — the markers ai and a2 can be put approximately equal. This means 
that their difference will be small: 



a2 — ai 



Besides, due to charge quasi-neutrality, Jacobians with respect to ai and a2 must be equal 
each other (here we put riio(ai) = ?T-2o(a2) = 1 without loss of generality): 

det||9x/9ai|| = det||(9x/9a2||. 

As a result, the infinitesimal vector d(a,t) relative to the argument a occurs divergence 
free: ddi/dai = 0. 

Then, summing (|7.12|) and ( [7.13|) and considering the limit e — > 0, we obtain the 



Weber-type transformation coinciding with (|6.4| ): 

u(a,t) = uo(a) + Va$ + Bo(a) x curU S 



(7.14) 



where vectors uo(a) and S are expressed through the Lagrangian invariants Pi(a) and 
P2(a) and displacement d between electron and ion by means of relations [pT| : 



u(a,t) = — (pi(a) + p2(a)), 

m 

d = curlaS. 

e 



Important that in ( [7. 14|) all terms are of the same order of magnitude (zero order relative 



to e). Curl of vectors pi(a) and P2(a2) yield the corresponding Cauchy invariants (|7.6|) 
and (^). 



8 Relabeling symmetry in MHD 

Now let us show how existence of new Lagrangian invariants corresponds to the relabeling 
symmetry. 

Consider the MHD Lagrangian H], 



8n 
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where we neglect by contribution from electric field in comparison with that from magnetic 
field. Here e{p) is specific internal energy. 

In terms of mapping x(a, t) the Lagrangian is rewritten as follows ||22|| : 



Here density and magnetic field are expressed by means of relations 

p = 1/Jx, h = (ho(a)Va)x/J,, 

and 

J^(a, t) = det||9x/9a|| 

is the Jacobian of mapping x = x(a,i(;) and initial density is put to equal 1. Notice, that 
variation of the action with by the Lagrangian ( ^.11 ) relative to x(a) gives the equation of 
motion ( |6.2| ) (or the equivalent equation for vector u (|6.3|)). 

Due to the presence of magnetic field in the Lagrangian ( ^.1|) , the relabeling symmetry, 
in comparison with ideal hydrodynamics, reduces. If the first two terms in ( p.l| ) are 
invariant with respect to all incompressible changes a — > a(b) with J\i, = 1, invariance 
of the last term, however, restricts the class of possible deformations up to the following 
class 

(ho(a)Va)b = ho(b). 

For infinitesimal transformations 

a ^ a + rg(a) 

where r is a (small) group parameter the vector g must satisfy two conditions: 

divag = 0, curl4g X ho] = 0. (8.2) 

The first condition is the same as for ideal hydrodynamics, the second one provides con- 
servation of magnetic field frozenness. 

The conservation laws generating by this symmetry, in accordance with Noether the- 
orem, can be obtained by standard scheme from the Lagrangian (|8.1|) . They are written 
through the infinitesimal deformation g(a) as integral over a: 

/ = y'(u,g(a))o!a (8.3) 

where the vector u is given by ( |2.6|) . Putting g = hg from this (infinite) family of integrals 
one gets the simplest one 




which represents a cross-helicity characterizing degree of mutual knottiness of vortex and 
magnetic lines. 

The conservation laws ( |8.3|) are compatible with the Weber-type transformation. Re- 
ally, substituting (|6^ ) into (|8.3|) and using ( ^.21) one leads to the relation 

J (uo(a),g(a))o?a. 
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Hence conservation of ( |8.3|) also follows. Note that if one would not suppose an indepen- 
dence of uq on t then, due to arbitrariness of g(a), this could be considered as independent 
verification of conservation of solenoidal field uq: 



dt 



uo = 0. 



The MHD equations expressed in terms of Lagrangian variables become Hamiltonian 
ones, as in usual mechanics, for momentum p = x and coordinate x. These variables 
assign the canonical Poisson structure. 

In the Eulerian representation the MHD equations can be written also in the Hamil- 



tonian form 14 



Pt = {p,H}, Vi = {v,iJ}, ht = {h,H}, 

where noncanonical Poisson bracket {F, G} is given by the expression (|1.9|) . As for ideal 
hydrodynamics, this Poisson bracket occurs to be degenerated. For example, the cross 
helicity Ich serves a Casimir for the bracket (|1.9|) . The reason of the Poisson bracket 
degeneracy is the same as for one-fluid hydrodynamics - it is connected with a relabeling 
symmetry of Lagrangian markers. 

For incompressible case the Poisson bracket ( |1.9| ) reduces so that it can be expressed 
only through magnetic field h and vorticity f2: 



curl-— X curl— — 
oh oil 



JG JF 
curl-— X curl— ■ 
oh oil 



dh 



(8.4) 



+ I \ n 

This bracket remains also degenerated. 



, , (5G 

cuhttt X curl-— 
oil oil 



dh. 



9 Variational principle for incompressible MHD 

By analogy with incompressible hydrodynamics, one can introduce magnetic line repre- 
sentation: 

h(r, t)= 5{y - R(a, t))(ho(a)Va)R(a, t)d^ai. (9.1) 

For vorticity the analog of vortex line parametrization (O) can be obtained, for instance, 
as a limit e ^ of the corresponding representations for the two- fluid system. Calculations 
give PT[ : 

0(r, t) = j 5{y - R(a, t))((flo(a) + curla[ho(a) x U(a, t)]) Va)R(a, t)d=^a, (9.2) 

Here the field U(a, t) is not assumed solenoidal, as well as the Jacobian of mapping 
r = R(a,t) is not equal to unity. 

From the corresponding limit of the two-fluid system to incompressible MHD it is 
possible also to get the expression for Lagrangian 



L = j t;3a([(hoVa)R X (UVa)R] ■ R0+ 



(9.3) 
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+1/3 J d^a{[Rt X R] • (f2oVa)R) -:h:{rj{R,U},h{R}}. 

The Hamiltonian of the incompressible MHD 'Hmhd in terms of U(a, t) and R(a,t) takes 
the form 

1 / ((ho(a)Va)R(a))2 



IMHD 

671 



- ~ ' det\\dR/da\\ ^ ^" 



1 /■ /■((ri(ai)Vi)R(ai)-(f2(a2)V2)R(a2)) ,3 ,3 



Svr 7 J |R(ai) - R(a2j 

where we introduce the notation 



c/^aid^a2, (9.4) 



n{a,t) = ^^o(a) + curla[ho(a) x U(a,t)]. 

Equations of motion for U and R follow from the variational principle for action with 
Lagrangian ( |9.3|) : 

[(hoVa)R X R,] ■ (dR/dax) = -5n/5Ux, (9.5) 

[(ri(a,t)Va)R X Rt] - [(hoVa)R x (UiVa)R] = 6n/6R. (9.6) 

These equations can be obtained also directly from the MHD system ( |1.5[ - |1.7|) by the same 
scheme as it was done for ideal hydrodynamics. 

Thus, we have variational principle for the MHD-type equations for two solenoidal 
vector fields. Their topological properties are fixed by rio(ci-) and ho (a). These quanti- 
ties represent Casimirs for the initial Poisson bracket (^.4|) . It is worth noting that the 
obtained equations of motion have the gauge invariant form. This gauge invariance is a 
remaining symmetry connected with relabeling of Lagrangian markers of magnetic lines 
in two-dimensional manifold which can be specified always locally. Coordinates of this 
manifold enumerate magnetic lines. This symmetry leads to conservation of volume of 
magnetic tubes including infinitesimally small magnetic tubes, namely, magnetic lines. 
This property explains why the Jacobian of the mapping r = R(a,t) can be not equal 
identically to unity. 
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